Enhanced vortex pinning by a composite antidot lattice in a superconducting Pb film 
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The use of artificial defects is known to enhance the superconducting critical parameters of thin 
films. In the case of conventional superconductors, regular arrays of submicron holes (antidots) 
substantially increase the critical temperature Tc{H) and critical current Ic{H) for all fields. Using 
electrical transport measurements, we study the effect of placing an additional small antidot in the 
unit cell of the array. This composite antidot lattice consists of two interpenetrating antidot square 
arrays with a different antidot size and the same lattice period. The smaller antidots are located 
exactly at the centers of the cells of the array of large antidots. We show that the composite antidot 
lattice can trap a higher number of flux quanta per unit cell inside the antidots, compared to a refer- 
ence antidot film without the additional small antidots in the center of the cells. As a consequence, 
the field range in which an enhanced critical current is observed is considerably expanded. Finally, 
the possible stable vortex lattice patterns at several matching fields are determined by molecular 
dynamics simulations. 



I. INTRODUCTION 

During the last decade, compelling evidence has shown 
that the introduction of an array of micro-holes (anti- 
dots) in a superconducting film has a profound influence 
on both the critical currenlji^ Ic{H) and the critical 
temperature^ Tc{H). Typically, at temperatures used 
for transport measurements, the antidots are able to trap 
only one flux quantum $o before saturation sets in. In 
this case, after the first matching field Hi , interstitial vor- 
tices appear in the sample, creating a "composite vortex 
lattice" , where part of the vortices is strongly pinned at 
the antidots and the rest occupies interstitial positions in 
between the antidots.-: Due to their higher mobility, the 
presence of interstitial vortices lowers the critical current 
Ic{H) significantly and broadens the R{T) transition. 

In this work, we study a composite antidot array, con- 
sisting of two interpenetrating square lattices with the 
same period d = 1.5 /im, but different antidot size 
(ai = 0.55/Ltm and a2 = 0.25 ^m). The two sublattices 
are shifted with respect to each other by half a unit cell 
along X- and y-directions, so that the small antidot is 
placed in the center of the unit cell of the lattice of large 
antidots (see Fig.nj. This arrangement of antidots cor- 
responds to the vortex lattice configuration at the second 
matching field in a sample with a single square array of 
antidots with rig = 1. The purpose of this experiment 
is to enlarge the field range where an enhanced criti- 
cal current Ic{H) is observed, by having efficient pin- 
ning sites exactly at the locations where the interstitial 
vortices would appear if the smaller antidots were not 
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present. 

This paper is organized as follows. In the next section 
we present some details concerning the sample prepa- 
ration and characterization. In section ITTll we study the 
phase boundary Tc {H) which allows us to identify the dif- 
ferent vortex configurations. Section llVl is devoted to the 
fiux pinning properties of this composite array of holes by 
measuring the critical current as a function of field and 
temperature. We show that a composite antidot lattice 
considerably increases the critical current at high fields. 
This effect results from the fact that for H > H2 the sat- 
urated small holes force the incoming vortices to occupy 
the big antidots. This situation persists until H = H4 
where the big antidots saturate and interstitial vortices 
form a more complex pattern. Finally, in Section we 
determine the most stable vortex patterns by means of 
molecular dynamics simulations. 



II. SAMPLE PREPARATION AND 
CHARACTERIZATION 



The used sample is a 50 nm-thick Pb film with a com- 
posite antidot lattice. The results obtained with this 
sample are directly contrasted with those measured on 
a reference antidot sample without the small holes, i.e. 
fli = 0.5 /J,m and 02 = ^m. In both cases, the bridge 
made for transport measurements (see Fig. ^a)) has a 
width of 300 /im and a voltage contact separation of 
2 mm. The unit cell of the composite antidot array is 
shown schematically in Fig. Q^b). The procedures fol- 
lowed to grow the samples are described in Ref. 5]. The 
transport measurements were carried out in a commer- 
cial PPMS-Quantum Design device with a temperature 
stability better than 0.5 mK. All measurements were per- 
formed with the field H applied perpendicular to the sur- 
face of the film. The critical temperature Tcq — 7.207 K 
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was determined from the resistive transition i?(T) in 
zero field, using a criterion of 10% of the normal state 
resistance i?„. 



FIG. 1: Layout of the Pb Him with a composite array of square 
antidots of two different sizes, (a) Geometry of the sample 
showing the patterned area in dark gray, (b) Schematic pre- 
sentation of a unit cell of the antidot array, (c) Atomic force 
micrograph of a 5 x 5 iiia^ area of the composite antidot ar- 
ray. The lattice period d is 1.5 /xm, the antidot sizes are 
ai = 0.55 fj,m and 02 — 0.25 pim. 



Due to the lateral nanostructuring, the effective width 
of the sample is reduced from 300 /im to 140 fim. Here, 
we have assumed that we can model the antidot sam- 
ple as a set of 200 parallel strips of width 0.7 /xm (= 
2 • 0.35 fj,m). This effective width was employed to cal- 
culate the resistivity p(7.5 K)= 5.33 • 10~* flm from the 
resistance i?(7.5 K)= 15.2 51. Using the listed valued 
for p£ = 4.88 • 10~^^J7m^, this resistivity value gives an 
elastic mean free path of ^ = 9 nm, and therefore a super- 
conducting coherence length ^(0) = 25 nm (in the dirty 
limit). These values are noticeably smaller than those 
obtained for the reference antidot sample {£ = 27 nm 
and ^(0) ~ 40 nm). Since in a film without antidots 
coevaporated with the sample containing the composite 
antidot was obtained £ — 27 nm, this difference seems to 
be caused by the more complex lift-off procedure due to 
the presence of the small holes. 

Knowing the mean free path £ and using the London 
penetration depth for the bulk^ Pb we obtain A(0) = 
71 nm. Due to the perforation, the effective penetra- 
tion depth increases, and therefore A should be modified 
according to^ 



A,(0) = , ^ ' - 86 nm , (1) 



where Sa and St are the area of the holes and the total 
area per unit cell, respectively. As a result, the Ginzburg- 
Landau parameter k amounts to k = = 3.4 > 
and therefore this sample is a Type-II superconductor.^*^ 

The sample has been characterized by means of atomic 
force microscopy. An AFM topograph of a 5 x 5 fim^ 
area of the film containing a composite antidot lattice is 
shown in Fig.^c). The root-mean-square roughness on 
a 1 fim.^ area of the sample in between the antidots is 
o'RMS = 3 nm. This value is about two times larger than 
for the plane film and the reference sample with antidots. 
This difference reinforces the idea that the film with the 
composite antidot lattice has suffered a small degradation 
due to a more complicated lift-off procedure. 



III. SUPERCONDUCTING Tc{H) PHASE 
BOUNDARY 

A. Experimental results 

We have measured the critical temperature Tc{H) as 
a function of field for the sample with a composite an- 
tidot lattice. The results obtained with a resistance cri- 
terion of 10 % of the normal state resistance i?„ and a 
measuring current of /ac — 10 fiA are shown in Fig. |21 
together with the phase boundary obtained for the ref- 
erence antidot film. The solid line depicts the expected 
upper critical field boundary of a plain film with the same 
coherence length as the reference antidot sample accord- 
ing to Hc2 = *i'o/2vr^(T)^. It is important to notice that 
the measured boundary of the reference antidot film is 
close to the Hcs (T) line corresponding to the surface nu- 
cleation of superconductivity around the holes, whereas 
the solid line represents the bulk superconducting tran- 
sition Hc2{T). As a rule, the experimentally determined 
critical temperature Tc{H) of a patterned sample turns 
out to be higher than that obtained for a plain film with 
the same coherence length^i 

FIG. 2: Tc{H) phase boundary for the film patterned with 
a composite antidot array, measured with an ac current of 
lac = 10 fiA and a resistance criterion of -Rcrit = 10% R„ 
(filled symbols). The open symbols show the phase bound- 
ary obtained for the reference antidot sample using the same 
criterion. The solid line is the calculated linear Tc{H) phase 
boundary for a plain film with the same coherence length 
^(0) — 40 nm as the antidot patterned film. The field axis is 
normalized to the first matching field Hi = 9.2 G. The tem- 
perature axis is normalized to Tco, the transition temperature 
at H = 0. 

Due to the presence of the antidot array, matching 
features appear in Tc{H) with a periodicity of Hi = 
$0/*^^ = 9.2 G, corresponding to the lattice parame- 
ter d — 1.5 /im. Although not all of them are very 
pronounced, local maxima are visible in the Tc{H) of 
the composite array for all integer matching fields Hn 
(n = 1,2, ...,6), whereas no evidence of rational matching 
features is observed. Thus, the addition of the extra an- 
tidot in the center of the unit cell of the array with large 
antidots, leaves the matching period unchanged. This is 
an important observation, since the composite antidot 
lattice can also be regarded as a square lattice, tilted by 
45°, with a unit cell twice as small as that of the original 
lattice. If this were the periodicity felt by the vortices, 
the matching period would amount to 18.4 G, which is 
twice as large as the observed period. In that case, one 
would expect the local maxima at even matching fields 
Hn (n — 2,4, ...) in Fig. |21to be more pronounced than 
the ones at odd matching fields Hn {n = 1,3, ...). Since 
this is not the case, we conclude that all these peaks cor- 
respond to integer matching fields, indicating that the 
main period felt by the vortices is the period of the lattice 
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with large antidots. 

In order to identify the vortex patterns at the match- 
ing fields, wc have plotted the R{T) transition width 
AT,{H)= r,(i?crit = 97%i?„) - Te(i?crit = 0.1%i?„) as 
a function of H in Fig. O (filled symbols). In this plot 
three different regimes can be clearly distinguished. For 
H < H4, the coherence length is larger than the width 
of the strands thus leading to a parabolic background in 
the Tc{H) phase boundary. In this so-called "collective" 
regime, we observe that the R{T) transition width re- 
mains almost constant. For fields higher than H^, an 
increase of the transition width can be observed, super- 
posed with matching features at H5 and Hq. We in- 
terpret the sudden increase in the transition width as a 
crossover to the regime where interstitial vortices appear 
in the sample. 

The interstitial regime is indicated by the gray area in 
Fig. 13 for the composite array. This regime ranges up 
to 3.6 ^(T) = d — a, i.e. up to ~ i/s, where a change 
in the ATc{H) slope can be observed. For higher fields, 
the single object regime is entered, where a linear phase 
boundary slightly distorted by an oscillation with period^ 
H* = ^o/oi ~ 69 G, is expected. Although the linear 
phase boundary is indeed observed, single object oscil- 
lations are difficult to resolve in the narrow field range 
investigated. For comparison, in the same figure we show 
Ar,(H)= Te(i?crit = 99%i?„) - T,(i?„it = 0.1%i?„) for 
the reference antidot sample (open symbols). From this 
curve we can infer that if the smaller additional antidots 
are absent, the crossover to the interstitial regime occurs 
at H ^ 1.5 Hi. Therefore the presence of the additional 
smaller antidots has substantially delayed the appearance 
of interstitial vortices. From the ATc{H) curve, we thus 
conclude that the total number of trapped flux quanta per 
unit cell of the antidot lattice is at least four. 



FIG. 3: Filled symbols: transition width A Tc{H)^ 
Tc(Rcrit = 97%R„) - Tc{RcTit = 0.1%R„) of the film with 
a composite antidot array, measured with a current of /ac ~ 
10 fiA. The gray box marks the "interstitial" regime, where 
the increase of ATc{H) indicates the appearance of inter- 
stitial vortices. Open symbols: transition width A Tc{H) — 
Tc{Rcrit = 99% Rn) -Tc(Rcrit = 0.1%i?„) of the reference an- 
tidot film for Jac = 10 fj.A. The thin black arrows indicate the 
ordinate scale for each curve. 



In other words, the observed saturation number is 
much larger than in the reference antidot sample, where 
only a <E>o-vortex can be pinned per antidot. Even taking 
into account the fact that the large antidots in the com- 
posite antidot lattice are slightly larger (a — 0.55 /im) 
than in the reference antidot sample (a — 0.5 fim) and 
the difference in coherence length, this is still a rather 
surprising observation. Indeed, the addition of one small 
hole per unit cell of the antidot array leads to an unex- 
pected increase of the number of pinned flux quanta per 
unit cell from one to four. 



B. Discussion 

To determine the number of flux quanta located in the 
large and the small antidots of the composite antidot ar- 
ray, a further investigation of the Tc{H) phase boundary 
is needed. As explained in Ref.|3], the background of the 
phase boundary is parabolic as long as 1.8 £,{T) > d — a 
(collective regime). In a square antidot lattice, this en- 
velope is described bjiii 



Tc{H) - T, 



cO 



\/3$o 



(2) 



By fitting the Tc{H) data points at the integer matching 
fields with this formula, the effective width d — a of the 
strands between the antidots can be deduced. 



FIG. 4: Tc{H) of the sample with a composite antidot ar- 
ray (open symbols). The solid line is a fit of the parabolic 
background (Eq. Q) shifted in field by one matching period. 
This line nicely interpolates between the Tc{Hn) points at the 
four matching fields (n = 2, 5). 

This procedure is, of course, not a priori valid for 
a composite antidot lattice. However, part of the 
Tc(H) phase boundary clearly shows a parabolic back- 
ground. The parabola providing the best envelope for the 
data points between H2 and H/^ is depicted by a solid line 
in Fig. 21 For fields higher than iJs, the measured curve 
deviates from the fitted parabola, as is expected for the 
"interstitial" regime starting for fields higher than H^. 
Strikingly, the fitted parabolic background turns out to 
be shifted by one matching period, having its maximum 
at Hi instead of at iJ = 0. By allowing this shift, the 
following parabolic envelope can be found for the second 
up to the fourth (or even fifth) matching peak, described 
by 



Tc{H) = 



cO 



1 - 



Mo(g-gi)^(0)7r(d-aeff) 



(3) 

using r^'g = 7.195 K and an effective antidot size of 
Ooff = 0.72 ^m. From this, we deduce that, after the 
second matching field H^ , the film with a composite anti- 
dot lattice with a filled small antidot behaves as if it would 
have contained a single square antidot array, but with a 
larger antidot size (ocff = 0.72 /im). 

The following scenario might explain such a behavior. 
Up to Hi , the vortices will be attracted towards the large 
antidots. Between Hi and H2, vortices begin to occupy 
the small antidots. Due to their size, these small antidots 
trap at most a single quantum vortex. They will there- 
fore be completely saturated at H2, creating a repulsive 
potential at the position of the small antidot. Fig. 
shows a schematic evolution of the potential landscape 
along a diagonal of the array (see dotted line in the in- 
set) that would be experienced by a vortex for H = 0, 
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ifi, and H2. Since the large antidots pin one flux quan- 
tum, at H — Hi a surface barrier has emerged at the 
antidot edges. For H — H2, the contribution to the po- 
tential of the small antidot at the center of the unit cell 
is strongly repulsive. When additional vortices enter the 
sample, they will be pushed towards the large antidots, 
leading to an increase of their effective saturation num- 
ber. In other words, the additional repulsive potential at 
the small antidots helps to increase the saturation number 
of the larger antidots. Within this scenario, for H > H2 , 
the phase boundary Tc{H) of the composite antidot lat- 
tice with a saturated smaller antidot resembles strongly 
the phase boundary expected for a simple square antidot 
lattice, without small antidots in the center, but with a 
larger antidot size. Due to this larger effective size, these 
antidots are then able to trap more vortices. We there- 
fore conclude that of the four flux quanta trapped per unit 
cell of the composite antidot lattice, one is pinned by the 
small antidot, while three are pushed into the larger holes. 
The net effect of the addition of the smaller hole in the 
antidot film, is to increase the effective pinning capacity 
(or effective saturation number) of the lattice with large 
antidots. This leads to a substantial broadening of the 
field range where a strong Tc{H) enhancement is ob- 
served. A similar picture was introduced by Doria and 
co-workers to explain the multiple trapping of vortices 
at high fields, as a result of the pressure exerted by the 
external vortices into the pinning site^ 

FIG. 5: Schematic representation of the potential along a di- 
agonal of the composite antidot array (see inset), experienced 
by a vortex entering the sample for H = 0, Hi, and H2. 

The potential at H = H2, drawn schematically in the 
lower panel of Fig. |31 can be seen as a checkerboard 
pattern with consecutively a repulsive and an attractive 
site. In a recent calculation, Lindquist and Riklund have 
modified the classical problem of a two-dimensional elec- 
tron gas exposed to a magnetic fieldM by adding a peri- 
odic checkerboard-like on-site potential with alternating 
signs^'^. Since the lowest energy level Ell{H), found for 
this classical problem, corresponds to the Tc{H) phase 
boundary of a simple square antidot lattice in the "collec- 
tive" regime, we believe that the addition of the checker- 
board on-site potential in the theory corresponds to our 
experimental situation for the composite antidot lattice 
at fields higher than the second matching field {H > H2). 
These authors showed that both the integer and the ra- 
tional matching features in Ell{H) are smeared out in 
the checkerboard system. These results may explain why 
the matching features in Tc{H) are much weaker for 
H > H2, than for H < H2 (see Fig.QJ. Indeed, it is only 
when H2 is exceeded that the checkerboard-like pinning 
potential is realized experimentally. 

Summarizing the results of the Tc{H) phase bound- 
ary measurements, the composite antidot lattice shows 
enhanced pinning, with many integer matching features 
appearing for fields up to Hq. The broadening of the 



R{T) transition after H4 indicates that at least four flux 
quanta can be trapped per antidot. From the shift in 
magnetic field of the parabolic background of Tc{H) with 
one matching period, we have deduced that the small an- 
tidot pins a <I>o-vortex, while 3$o-vortices are trapped by 
the large antidots. The presence of the additional small 
antidot in the center of the unit cell has therefore led 
to a substantial broadening of the collective regime, or, 
in other words, to an increase of the effective saturation 
number Us of the large antidots from one to three. 



IV. CRITICAL CURRENT AS A FUNCTION OF 
FIELD 

So far we have explored the normal-superconducting 
boundary in order to experimentally determine the exten- 
sion of the different regimes given by the ratio ^(T)/((i — 
a). Now we turn to isothermal critical current measure- 
ments which allow us to study the vortex dynamics deep 
in the superconducting state. 

The critical current as a function of field Ic{H) was 
measured using a 100 /iV voltage criterion for several 
temperatures close to Tc{H). The results, in order of 
decreasing temperature, are shown in Figs. El and [7| 

The absolute value of the critical current density at 
zero field for the composite antidot array amounts to 
Ico - Ic{H - 0) - 6.8 • 10^^ at T/Teo = 0.974. 
This value is a factor ^ 3 lower than the critical current 
density obtained for the reference antidot lattice, in part 
due to the difference in A(0) and ^(0), which eventually 
determine the pinning properties. The specific geome- 
try of the lateral nanopatterning might also influence the 
current distribution throughout the film, hereby also af- 
fecting the critical current Ico- In order to compare the 
pinning properties of the film with the composite antidot 
lattice (open symbols) and the reference antidot lattice 
(solid lines) measured at the same reduced temperature 
we have normalized the critical current by Ico- Notice 
that since the saturation number is mainly determined 
by the coherence lengthi^ii^ ^ (T) which in turn depends 
solely on the reduced temperature t = T/Tc, regardless 
the value of the applied magnetic field, it is enough to 
compare the results obtained on these samples at the 
same reduced temperature, without the necessity of nor- 
malizing the field. 

FIG. 6; Normalized critical current at T/Tco = 0.997 and 
T /Tco = 0.993 of a film with a composite antidot array. The 
curves were measured for H < (open symbols) and sym- 
metrized for clarity for H > (dashed line). For comparison, 
the solid line shows the normalized critical current Ic{H)/Ico 
obtained for the film with the reference antidot lattice. 

The Ic{H)/Ico curves for the film with a composite an- 
tidot array (Figs. Eland 01 have been measured for nega- 
tive fields (open symbols) and symmetrized for H > for 
clarity (dashed lines). All curves show distinct periodic 
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FIG. 7: Normalized critical current at T/T^o = 0.986 and 
T /Tco = 0.974 of a film with a composite antidot array. The 
curves were measured for H < (open symbols) and sym- 
metrized for clarity for H > (dashed line). For comparison, 
the solid line shows the normalized critical current Ic(H)/IcO 
obtained for the film with the reference antidot lattice. 



matching features, with a period Hi = ^ ~ 9.2 G corre- 
sponding to the unit cell of the lattice with the large (or 
the small) antidots (d = 1.5 fim). As we pointed out be- 
fore, the periodicity felt by the vortices is defined by one 
of the interpenetrating antidot lattices of the composite 
array, and not by the resulting square lattice with a unit 
cell which is twice as small and rotated by 45°. 

In the upper panel of Fig. {T/Tco = 0.997), the 
Ic{H)/Ico curve of the film with a composite antidot lat- 
tice shows a peak-like behavior with sharp maxima at 
Hi, H2, and H3. This behavior is expected at tempera- 
tures sufficiently close to Tco, where it is not possible to 
have interstitial vortices in the superconducting strands 
between the antidots. As we have already deduced from 
the shape of the Tc{H) phase boundary (Section IIII|I . 
interstitial vortices indeed appear in the sample only for 
T/Tco < 0.994. 

At a lower temperature, T/Tco = 0.993 (Fig. El lower 
panel), and all temperatures below that (Fig.EJ, a strong 
enhancement of Ic{H)/Ico in the film with a composite 
antidot lattice can be found for fields higher than the 
first matching field Hi, compared to the reference anti- 
dot lattice. The reason for this lies in the ability of the 
composite antidot lattice to pin more flux quanta inside 
the antidots compared to the reference antidot array (see 
Section Iml. 

It should be noted that the field range where the film 
has a finite critical current, i.e. where the film remains 
superconducting, is considerably broader for the compos- 
ite than for the reference antidot array. 

The appearance and sharpness of the matching fea- 
tures in the Ic{H)/Ico curves, are temperature depen- 
dent. At T/Tco = 0.993 (Fig. lower panel), every in- 
teger matching peak up to Hq can be clearly seen. The 
maxima at Hi, H2, and H3 are very pronounced. At 
H4 and iJs, one finds cusps rather than local maxima in 
Ic{H)/ Ico- The matching feature at Hq is again peak- 
like. This indicates that the vortex patterns formed at 
H4 and Hq are less stable than the vortex configuration 
at Hq. 

When the temperature is lowered to T/Tco = 0.986 
(Fig. upper panel), we find again sharp matching fea- 
tures in Ic{H)/Ico at Hi, H2 and H^, and only very weak 
cusps at Hi and H^. At Hq, the local maximum has de- 
veloped into a pronounced cusp, after which a substantial 
change in the Ic{H)/Ico slope occurs. A second smaller 
slope change can be found at H-;. At the lowest mea- 
sured temperature, T/T^q = 0.974 (Fig.Q lower panel), 
the only matching features left are the sharp local max- 
ima at Hi , H2 , and H^ , and one pronounced cusp at iJy . 



FIG. 8: Suggested vortex pattern at Hi, H5, Hq, and Hj. 
All patterns have been obtained by molecular dynamics sim- 
ulations by an annealing procedure, except the one at H4,. 
Open circles and black dots represent pinning sites and single 
quantum vortices, respectively. 

It appears that, at this temperature, the seventh match- 
ing field Hj plays the same role as the sixth matching 
field Hq at T/Tco = 0.986. This fact leads us to believe 
that at T/Tco = 0.974, the total number of trapped flux 
quanta per unit cell of the composite lattice, has increased 
from four to five. 

It is worth to notice that the normalized critical cur- 
rent at the first matching field Hi reaches approximately 
the same value for the film with the composite and with 
the reference antidot lattice, except for the Ic{H)/Ico 
curve taken at T/Tco = 0.997. This fact makes the film 
with the reference antidot array a good candidate to com- 
pare its pinning properties with those of the composite 
antidot array. 



V. STABLE VORTEX PATTERNS 

The periodic matching features in the Tc{H) phase 
boundary and in the critical current versus field curves 
Ic{H) demonstrate that the composite antidot lattice 
can stabilize commensurate vortex lattices at several 
magnetic fields. From the results presented in Section lTTTI 
it is clear that the large antidots trap at least three flux 
quanta, while the smaller antidots are saturated after 
pinning one single quantum vortex. The vortex patterns 
expected in the composite antidot lattice will therefore 
differ from the known patterns in simple square pinning 
arrays. 

We have performed molecular dynamics simulations to 
obtain the vortex patterns at the matching fields H^,, Hq, 
and Hj. To model the composite vortex lattice, two in- 
terpenetrating arrays of Gaussian sites with a different 
radius and a different pinning force were used. This was 
necessary to obtain an occupation of 3$o-vortices in the 
large pinning sites, and of one $o-vortex in the small 
pinning sites. By applying an annealing course, in which 
the temperature is lowered, starting from a high tem- 
perature and a random distribution of vortices, the most 
stable configuration in the given pinning potential can be 
found. As the temperature is lowered, the pinning sites 
become smaller and stronger, scaling with ^(T) in the 
following wayi^: 



Fp cx Fpocxp (^-^j , (4) 

where FpQ is the single site pinning force. Additionally, 
the vortex- vortex interaction range reduces with decreas- 
ing temperature, due to the decrease of the magnetic 
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penetration depth A. In this type of simulation, the oc- 
cupation of the pinning sites lowers with decreasing tem- 
perature. The annealing course was therefore stopped 
when the occupation, corresponding to our experimen- 
tal situation, was achieved. Figure |21 shows the vortex 
configurations we suggest for H4, H5, Hq, and H7. All 
patterns, except the one at H^, are obtained from the 
molecular dynamics simulations. Multi-quanta vortices 
are represented in this model by a multiple occupation 
of a pinning site with (repulsive) single-quantum vortices. 
Since in the experiment, the pinning sites consist of real 
holes in the film, the vortices trapped in the same pinning 
site will be interpreted as multi-quanta vortices, even 
though they are depicted as separate single flux quan- 
tum entities in the plots. Actually, this model can be 
experimentally realised by an array of non fully perfo- 
rated (or blind) holes. Reported results on such systems 
showed that blind holes are weaker pinning centers than 
antidots, although the overall features in both cases are 
very much alike ^ 

The vortex pattern at H4, which is drawn schemati- 
cally and was not calculated, depicts all antidots occupied 
with the maximum number of vortices. The large anti- 
dots trap 3<i>o-vortices, the smaller antidots trap a $0- 
vortex. No interstitial vortices are present in the sample. 
In this case, one would expect the matching feature at 
H4 to be of the same kind as the one at H3. Surprisingly, 
the Ic{H) curves (see e.g. Fig.[71l show only very weak 
matching features at H4. However, the Tc{H) measure- 
ments and the fact that the matching peak in Ic{H) at 
Hq is very well defined (see discussion below), both indi- 
cate that there should be four vortices trapped per unit 
cell of the antidot array, leading to the suggested vortex 
pattern. 

At i/s , there is one interstitial vortex present per unit 
cell of the array. It is, however, not evident where this 
vortex is located, since the most logical position, at the 
center of the unit cell, is already occupied by the filled 
smaller antidot. One can see a tendency of the interstitial 
vortices to form diagonal lines, which make zigzag traces 
across the sample, indicated by dashed lines in Fig.|Hfb). 
However, the long range order which is needed to make a 
regular pattern, with for example straight diagonal lines 
or a regular zigzag, is lacking at this field. The pattern 
found in molecular dynamics simulations for H^, is con- 
sequently not very stable. 

At Hq, a highly symmetric vortex pattern is formed. In 
this case, two interstitial vortices are present per unit cell, 
which are positioned approximately at the center of the 
line connecting two neighboring large antidots. Due to its 
high symmetry, the vortex pattern at ffg is very stable. 
Remarkably, the calculations for i/gj resulting in a very 
regular pattern, have been performed under the same 
conditions as the ones at iJs, where no regular pattern 
could be found. This is an indication of the different 
stability of the vortex patterns at and iJg- Indeed, 
the matching feature at Hq in Tc{H) or Ic{H)/Ico is 
always more pronounced than at (see e.g. Figs. ^ 



andUl upper panel). For the Ic{H)/Ico curve measured 
at T/Tco = 0.974 (Fig. □ lower curve), the matching 
cusp at H'! becomes rather sharp. We suggest that at 
this temperature, the large antidots are able to trap four 
flux quanta instead of three. In that case, the expected 
vortex pattern for i/y (see Fig.|Sl), resembles the pattern 
calculated for the sixth integer matching fleld, but with 
four flux quanta occupying the large antidots instead of 
three. 

To obtain a regular pattern at i?7, with three flux 
quanta pinned in the large antidots, three interstitial vor- 
tices have to be accommodated per unit cell of the com- 
posite array. The calculations were not able to produce a 
regular vortex pattern with an occupation of three at the 
large pinning sites and one at the smaller pinning sites. 
This indicates that the stability of a vortex pattern at 
i?7 is not very high. On the other hand, the calculation 
method, where 3 ^o-vortices are represented by three 
separate <l>o-vortices, might in this case also affect the 
outcome of the simulation. Indeed, the cylindrical sym- 
metry of the 3 $o-vortex, as it occurs in the experiment, 
will make it easier to obtain a regular vortex pattern. 

The vortex patterns suggested in this section remain 
to be directly verified by using a local scanning tech- 
nique, such as low temperature scanning Hall probe mi- 
croscopy. Further insight into the vortex pinning and dy- 
namics in systems with a composite pinning array may 
also be gained from molecular dynamics simulations. 



VI. CONCLUSIONS 

We have used a composite antidot lattice, consisting of 
two interpenetrating antidot arrays with a different an- 
tidot size, but with the same lattice period, as a strong 
periodic pinning potential for the vortex lattice in a su- 
perconducting film. The shift between the two lattices 
is such that the smaller antidots are situated exactly at 
the centers of the cells of the array of large antidots. We 
have shown that this pinning array can stabilize the vor- 
tex lattice at several matching fields from Hi to iJy. 

Measurements of the critical temperature Tc{H) and 
current Ic{H)/Ico as a function of magnetic field, have 
demonstrated that the composite antidot lattice can trap 
a considerably higher amount of flux quanta per unit cell 
(four or flve instead of one) inside the antidots, com- 
pared to a reference antidot fllm without the additional 
small antidots in the center of the unit cell. This means 
that the appearance of interstitial vortices in the compos- 
ite antidot lattice is delayed to higher magnetic flelds. 
The presence of the smaller antidots has therefore in- 
creased the effective saturation number of the large an- 
tidots, which has led to a considerable expansion of the 
fleld range in which an enhanced critical current is ob- 
served. 
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